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Abstract: The quantization of mirror curves to toric Calabi-Yau threefolds leads to trace class 
operators, and it has been conjectured that the spectral properties of these operators provide a 
non-perturbative realization of topological string theory on these backgrounds. In this paper, we 
find an explicit form for the integral kernel of the trace class operator in the case of local P 1 x P 1 , 
in terms of Faddeev’s quantum dilogarithm. The matrix model associated to this integral kernel 
is an 0(2) model, which generalizes the ABJ(M) matrix model. We find its exact planar limit, 
and we provide detailed evidence that its l/N expansion captures the all genus topological string 
free energy on local P 1 x P 1 . 
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1 Introduction 


Topological strings on Calabi-Yau (CY) manifolds, just like other string theories, are only defined 
in perturbation theory, in terms of a genus expansion. In the closed string sector, the topological 
string free energies compute the Gromov-Witten invariants of the CY target, and for this reason 
topological string theory has played a prominent role in the interface of string theory, geometry, 
and mathematical physics. 

Recently, it has been conjectured in [1] that topological strings on toric CY threefolds are 
captured, non-perturbatively, by the spectral theory of quantum-mechanical, trace class opera¬ 
tors. These operators arise naturally in the quantization of their mirror curves. The conjecture 
of [1] builds upon previous ideas on quantization and mirror symmetry [2-5], but it also incor¬ 
porates many conceptual aspects of large N dualities. In fact, many of the crucial ingredients in 
the proposal of [1] were first unveiled in the study of the ABJM matrix model at large N [6-11] 1 . 
As spelled out in detail in [13], one way of formulating the spectral theory/mirror symmetry 
correspondence of [1] is by considering the so-called fermionic traces Z(N,h ) of the trace class 
operator (see section 2.1 for a precise definition.) It turns out that, in the ’t Hooft limit, 


N 


00 , 


h —> 00 , — = A fixed, 

n 


( 1 . 1 ) 


1 The approach of [10] was first applied to topological string theory in [12], but it requires corrections which are 
incorporated in [1[. 
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these traces have an asymptotic expansion of the form, 

logZ(N,H) = J2j 7 g(\)fr 2 - 29 - (1.2) 

g> o 

According to the conjecture of [1], the functions F g { A) should be the genus g free energies of the 
standard topological string, in the so-called conifold frame. The ’t Hooft parameter A is a flat 
coordinate for the CY moduli space and is given by the vanishing period at the conifold point 
(we are assuming here that the mirror curve has genus one.) In this way, the weakly coupled 
topological string emerges in a limit in which the quantum-mechanical problem is strongly coupled 
(since h —> oo). On the other hand, the double-scaling limit 

N —> oo, h -» 0, Nh = g fixed, (1-3) 

corresponds to the WKB expansion in the quantum-mechanical problem, and it is captured by the 
Nekrasov-Shatashvili (NS) limit of the refined topological string, in agreement with the results 
of [4, 5]. An obvious corollary of (1.2) is that the fermionic spectral traces Z(N,h) of the trace 
class operator provide a non-perturbative definition of the topological string partition function, in 
the spirit of large N dualities. From a more physical point of view, one can regard Z(N,h) as 
the canonical partition function of a quantum ideal gas of N fermions, where the operator plays 
the role of density matrix [6]. 

As explained in [13], if the kernel of the operator arising in the quantization of the mirror 
curve is known explicitly, then the fermionic spectral traces can be computed by a matrix model. 
Fortunately, it was shown in [14] that, for some simple mirror curves (leading to so-called three- 
term operators), one can compute the corresponding kernels in closed form, in terms of Faddeev’s 
quantum dilogarithm. This made it possible to verify the trace class property conjectured in [1]. 
Armed with these kernels, one can compute the fermionic spectral traces Z(N,h), which are 
given by a generalized 0(2) matrix model of the type considered in [15]. In [13] this matrix 
model was studied in the 1/N expansion, and it was checked in detail that, for local P 2 and a 
certain limit of local P 2 , (1-2) gives indeed the topological string free energies. 

This paper extends the results of [13, 14] to an important local CY, namely local P 1 x P 1 . 
In this case, the mirror curve has genus one, therefore one modulus, but it also involves a 
mass parameter, since the geometry has two Kahler parameters. The quantization of this curve 
leads to a four-term operator. By using the quantum pentagon identity for Faddeev’s quantum 
dilogarithm, we find an explicit expression for the integral kernel of the corresponding trace class 
operator. The matrix model obtained from this kernel turns out to be an 0(2) matrix model. 
We compute some spectral traces at finite N, as a function of the mass parameter, which agree 
with the predictions of the conjecture in [1], as shown in [16]. We also study the matrix model 
in the large N limit. This can be done by doing perturbation theory in the ’t Hooft coupling, 
as in [13], but we can also use the general techniques of [17, 18], as developed in [19], to solve 
exactly for its planar limit. We compare the resulting 1/N expansion with the topological string 
genus expansion, and we find a detailed agreement. 

It is known, geometrically, that the topological string on P 1 x P 1 is equivalent to the topo¬ 
logical string on local F 2 by a simple change of parameters [16] (this leads to a relation between 
the Gromov-Witten invariants of the two geometries, as pointed out in [16, 20]). We show that 
this equivalence holds as a unitary equivalence between the corresponding trace class operators. 
This allows us to extend all of our results to local P 2 with arbitrary moduli, extending in this 
way the analysis presented in [13]. 
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The trace class operator obtained by quantization of the mirror curve of local P 1 x P 1 can be 
regarded as a generalization of the density matrix appearing in the Fermi gas formulation of the 
ABJ(M) matrix model [6, 21-23], after an analytic continuation to complex mass parameters. 
Therefore, one can rederive from our results various aspects of the ABJ(M) matrix model. For 
example, we show that the exact planar solution of the 0(2) matrix model reproduces the planar 
free energy of the ABJ(M) matrix model obtained in [24]. 

It is natural to ask how our matrix model for local P 1 X P 1 compares to a previous proposal 
in [25]. This proposal is based on a generalization of the Gopakumar-Vafa large N duality [26], 
in which topological string theory on local P 1 X P 1 is described by large N, U(N ) Chern-Simons 
theory on the lens space L{2, 1) [25]. When this is combined with the results of [27], one obtains 
a matrix model description of topological strings on local P 1 x P 1 which has been studied in 
some detail [25, 28, 29]. There are however many important differences between these matrix 
models. First of all, the matrix model of [25] is a two-cut matrix model, while our model is a 
one-cut matrix model. This leads to important differences at the non-perturbative level, since 
in the model of [25] the two Kahler parameters of local P 1 x P 1 are discretized (they correspond 
to the two filling fractions of the two-cut matrix model), while in the matrix model described 
here only the “diagonal” Kahler parameter is discretized. Another difference between these two 
matrix models is that the weak’t Hooft coupling expansion of the model in [25] corresponds to 
the so-called orbifold point in the moduli space of local P 1 x P 1 , while in the model considered 
here it corresponds to the conifold point. Both points lead to logarithmic periods (which are 
in fact needed to match the Gaussian behavior of the matrix models), but they are different. 
It would be interesting to understand in more detail the relationship between the two matrix 
models, specially at the non-perturbative level, but we will not pursue this problem here 2 . Note 
that, if the conjecture of [1] is true, the matrix model description in terms of kernels of trace 
class operators studied in this paper is likely to apply to all toric CY threefolds. In contrast, the 
large N duality of [25] applies only to a special type of geometries, obtained as ADE quotients 
of the resolved conifold. 

This paper is organized as follows. In section 2 we elaborate on [14] and obtain an explicit 
representation for the integral kernel of the trace class operator associated to local P 1 x P 1 . We 
also write down an 0(2) matrix model computing the fermionic spectral traces, and we study its 
1 /N expansion. We obtain perturbative results as well as a closed form expression for the planar 
free energy, which can be expanded at both weak and strong coupling. In addition, we show 
how many known results for the ABJ(M) matrix model can be recovered from this solution. In 
section 3 we compare successfully the 1/N expansion of the matrix model with the topological 
string free energies of local P 1 x P 1 , which we compute around a generic point in the conifold 
locus. We conclude in section 4 and we list some open problems for the future. In the Appendix, 
we list some properties of the quantum dilogarithm which are used in section 2. 


2 Operators, kernels and matrix models 

2.1 Integral kernel and matrix model for local P 1 x P 1 

As explained in [1, 14], given the mirror curve to a toric CY threefold, one can quantize it to 
obtain a trace class operator p. Although this procedure can be followed for any toric geometry, 
the simpler case is that of toric (almost) del Pezzo CY threefolds, defined as the total space of 

2 We would like to thank R. Schiappa for raising this issue. 
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the canonical line bundle on a toric (almost) del Pezzo surface S, 


X = 0(K S ) (2.1) 

In this case, the mirror curve has genus one. The complex moduli of the curve involve a “true” 
geometric modulus u as well as a set of “mass” parameters nii, i = 1, • • • , r, where r depends on 
the geometry under consideration [30, 31]. The mirror curves can be put in the “canonical” form 

W(e x ,e y ) = O s (x,y) + u = 0, (2.2) 


where Os(x,y) is given by 


k+2 

Os(x,y) = ^2 ex P + u 2 )r y + > ( 2 -3) 

i=l 

and fi(nij) are suitable functions of the parameters nij. The vectors \ can be obtained from 
the toric description of the CY threefold. The mirror curve (2.2) is quantized by standard Weyl 
quantization. In particular, x. y are promoted to self-adjoint Heisenberg operators x, y, satisfying 
the commutation relation 

[x,y] = i H, (2.4) 

and ordering ambiguities are resolved by Weyl’s prescription. In this way, Os(x,y ) becomes an 
operator, which will be denoted by 0 5 . As conjectured in [ 1 ] and proved in [14] in many cases, 
the inverse operator 

PS = Og 1 (2.5) 

is of trace class. 

In this paper we will focus on the important local del Pezzo CY threefold in which S = 
P 1 x P 1 = Fo, and usually called local P 1 x P 1 or local Fo- Topological string theory on this 
background is known to have various applications: it engineers geometrically SU( 2 ) Seiberg- 
Witten theory [32], it is dual to Chern-Simons theory on the lens space L( 2,1) [25], and it is 
closely related to the partition function of ABJ(M) theory on the three-sphere [ 6 , 24]. In this 
case, the function Os{x,y) is given by, 

Of 0 (x, y) = e x + m Wo e~ x + e y + e~ y , (2.6) 

and depends on a mass parameter that we denote by mp 0 . In principle, we will take m^ 0 to be 
real and positive, but as we will see it is possible to extend some of the results to complex values 
of mw 0 . 

We would like to find an explicit expression for the kernel of the operator p^ 0 . As for the 
three-term operators analyzed in [14], this kernel will involve in a crucial way Faddeev’s quantum 
dilogarithm <I>b(®) [33-35], see the Appendix for its definition and some of its basic properties. 
In addition, the function < hb(-c) has the following features. If p and q are self-adjoint Heisenberg 
operators satisfying, 

[p,q] = (27ri)-\ (2.7) 

the quantum dilogarithm satisfies [14] 

$ b (p)e 2wbq 4>b(p) = e 27rbq + e 2 ^ b(p+q) , 

( 2 . 8 ) 

$£(q) $ b (p)e 2,rbq $b(p) $b(q) = e 27rbq + e 2 ^ b(p+q) + e 2wb(p+2q) . 
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One also has the important quantum pentagon identity [36] 


$b(p) $b(q) = $b(q) $b(p + q) $b(p)- 

The quantization of (2.6) leads to the operator, 

Of 0 = e x + mi 0 e“ x + e y + e _y . 


Let us set 

and 

By using (2.8), we find 


h = 7rb 2 


x = 7rb(p + 2q), y = 7rbp. 


(2.9) 


( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


e x/2 0 Fo e x/2 - m Fo = e 2x + e x+y + e x " y = e 2 ^ b(p+2q) + e 2?rb(p+q) + e 2 ^ bq 

= <T b (q) <h b (p)e 2wbq <T b (p) 4> b (q). (2.13) 

Therefore, 


<T b (p) d> b (q)e x / 2 0 Fo e x / 2 <T b (q) 3> b (p) = m Fo + e 27rbq = m Fo (l + e 2 ^"^) 

d>b(q ~ bg/7T - ib/2) 

d>b(q - b£/7T + ib/2) 


= niw 0 


, (2.14) 


where the parameter £ is related to ?n Fo through the equation 


m ¥o = e 2b ' 5 . 


(2.15) 


Let us now define the operator 

B = <h b (q - b£/ 7 r - ib/2) $ b (p) < h b (q)e 7 rbp/, 2 e 7rbq . (2.16) 

We obtain the following formula [14] 

°d = (2-17) 

On the other hand, we can use the quantum pentagon relation (2.9) to write the operator B as 

e —Ob/2)V(4^i) $b(p)B = $ * (p + q _ b?/7r _ ib/2) $* (q _ b?/?r _ jb/2) £ b(q)e7 rb(p+q)/2 e ^bq/2 

= $£(p + q - b£/ 7 T - ib/2)e " b(p+q)/2 $£(q - b£/vr - ib/4) $ b (q + ib/4)e 7rbq / 2 . (2.18) 
If we introduce new momentum and position operators by 

P 7 = P + q “ b£/vr, q' = q-b£/2vr, (2.19) 


1 

2cosh(7rbp / ) 


r(q'), 


we find 


P¥o = e b2?/2 /(q') 


( 2 . 20 ) 




where 


( 2 . 21 ) 


f(n) = r ^bg/ 2 $ b(q ~ b^/27r + ib/4) 

<h b (q + b£/27r — ib/4)' 

In the position representation for the operators p 7 and q 7 , we obtain the integral kernel, 

= <z.|0 F » = 


( 2 . 22 ) 


As shown already in [14], this is a positive-definite, trace class operator on L 2 (M). Note that 
Pf 0 (* t 1; x 2 ) is related by a unitary transformation to the symmetric, real kernel 


, —b 2 £/2 l/(ai)ll/(32)| 

2 bcosh (tt X i ^ X2 ) ’ 


(2.23) 


which is of the type considered in [37, 38]. In particular, as shown in these references, its diagonal 
resolvent can be obtained from a TBA-like system of non-linear integral equations. 

The spectral information of a trace class operator p depending on a parameter h and acting 
on a Hilbert space T~L can be encoded in different ways. The spectral traces of p are defined by 


Zi = Tr w p e , l = 1, 2, ■ ■ ■ 


(2.24) 


The fermionic spectral traces are given by 


Z(N, h) = Tv aN(H) (A n (p)) , N = 1, 2, • • • , (2.25) 

where the operator A N (p) is defined by p® N acting on A N (FL). The generating function of the 
fermionic spectral traces is the Fredholm or spectral determinant of p: 

OO 

S(k, K) = det(l + up) = 1 + Z{N, K)k n , (2.26) 

N=l 

and is an entire function of k due to the trace class property of p [39]. A well-known theorem 
of Fredholm (see chapter 3 of [39] for a proof) states that Z(N,h) has the matrix-model-like 
representation 

Z(N,h) = J d N x det (p(xi, Xj)) . (2.27) 

In this equation, p(x 1 ,^ 2 ) is the integral kernel of the operator p , 

p{xi,x 2 ) = {xi\p\x 2 ). (2.28) 

The spectral traces (2.24) and the fermionic spectral traces are closely related, since one has that 

J(k) = log a(K,h) = ~J (“«)* ■ ( 2 - 29 ) 

l= 1 

The above quantities can be interpreted, more physically, in terms of an ideal Fermi gas of N 
particles, as in [ 6 ]. In this setting, p is the canonical density matrix, Z(N,h) is the canonical 
partition function of the gas, E(n,h) is the grand canonical partition function, and J{k) is the 
grand potential. 
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Since we have an explicit formula for the integral kernel of p ¥o , we can write down an explicit 
expression for the integral (2.27). By using Cauchy’s identity, as in [6, 13, 40], 


rii<j 

2 sinh 


2 sinh 

(?) 


n,j 2 cosh 1 

j \ 



- = det 


ij 


fli Vj 
2 


2 cosh 

E(-D' w n 






ceSjv * 2 cosh 

we obtain the following matrix model representation for the fermionic traces of p Fo , 

n , 9 rr a oini, 2 ( Ui ~ u J 


Zw 0 (N, h) = 


5 -b 2 £7V/2 r d N 


U 


N\ J (2vr) 


N 


n|/(£) 

i= 1 v ' 


2 n ; <i 4sinh 


n i,j 2 cosh 


where the variables rq are related to the original variables Xi by 

2 7T 

Ui = --Xi. 

b 


(2.30) 


(2.31) 


(2.32) 


2.2 Relation to local F 2 and spectral traces 

It is known that topological string theory on the local F 2 geometry is closely related to topological 
string theory on local Fo [16]. It turns out that this equivalence also holds at the level of the 
corresponding quantum operators. To see this, let us first redefine the operators appearing in 
(2.10) as, 

h = 2irb 2 , x = 27rbq, y = 27rbp. (2.33) 

We then have, 


/-x'i pX/2 = ^bq Mp — q~ ib/2) 
4>b (P - q + ib/2) 


+ e y = e x/2 (1 + e y_x ) e x/2 = —-—^4^e 7rbq 


= ^b (p — q) e 27rbq 4> b (p — q)” 1 . (2.34) 


Therefore, 


$b(p-q) 


Of 0 <Mp - q) - e 2?rbq = - \ (m Fo e- 2 - bq + e" 2 ^) 4> b (p - q) 


>(p-q) 

= m Fo e~ 7rbq ^ b j P ~ q ~ l b ^j e -7rbq + e ~^ b P ^ b(p - q ~ l^ e"^ 


4>b(p - q + ib/2) 4> b (p - q + ib/2) 

= m Fo (e" 27rbq + e 2wb (P- 2q )) + e~ 2nbp + e~ 2nbq = (1 + m ¥o )e~ 2nbp + m ¥o e 2nb ( p ~ 2q ) + < 

or in terms of original variables 

0 Fo 4> b (p - q) = e x + (1 + m Fo )e _x + m Fo e y_2x + e _y . 


4 — 27 rbp 


^b(p - q) 


(2.35) 


(2.36) 
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By defining new variables 

X =x + u. y'= y — 2x-3 u, v = log(mF 0 ), (2.37) 

we rewrite (2.36) as follows 

^ b (p L _q) m F 0 1/lQF o $ b(p - q) = e x ' + (m^ 2 + mp o 1/2 )e“ x ' + e y ' + e“ 2x '“ y '. (2.38) 

We conclude that the operator m Fo 1//4 OF 0 is unitarily equivalent to the operator 

0f 2 = e x + m.F 2 e _x + e y + e~ 2x_y , (2.39) 

corresponding to the local F 2 geometry [1, 14], after the substitution 

?b F2 = ruj 2 + m~l /2 . (2.40) 

In the CY geometries, the rescaling by m^ 1 ^ 4 leads, in view of (2.2), to the following relation 
between the moduli, 

% 2 = ^ F 0 1/4 %0 • ( 2 ‘ 41 ) 

The relationships (2.40), (2.41) agree precisely with those found by a direct analysis of the 
topological string in these geometries [16]. This means in particular that any test of the conjecture 
of [1] for local Fo leads automatically to a corresponding test for local F 2 . The unitary equivalence 
of the two operators also leads to the following equality of spectral traces, 

Tr P¥ 2 ( m F 2 ) = ^o^PFo (™f 0 ), (2-42) 

after the substitution (2.40). 

Using the expression for the integral kernel in (2.22), as well as (2.42), we can in principle 
compute explicitly the first spectral traces. According to [1], we should expect simplifications in 
the so-called maximally supersymmetric case h = 2 -k, which corresponds to 

b = V2 (2.43) 


in (2.11). For this value of b, we can use the functional equation (A.9b) satisfied by the quantum 
dilogarithm to obtain the following expression in terms of elementary functions, 

1 


I fix ')I —- — 

4 COSh ( W2(x-bg/27T) ) cogh ^V2(x+b{/2n) 

After an appropriate change of variables, we obtain, 

T = —1/4 [°° _d u _ = 1 log(m Fo ) 

PF ° 87 t f ° J_ 00 cosh(u) cosh(rt - \/2b£/2) 87 r m V 2 _ ^ 


(2.44) 


(2.45) 


The second trace is a little bit more complicated. We find 


1 POO POO 


dudv 


-00 cosh(ti) cosh(it — y/2b^/2n) cosh(u) cosh(u — \Z2b^/2n) cosh(u — v ) 2 


m 


-1/2 

'Fo 


167T 2 


log(m Fo ) 


1/2 


, m F 0 m F 0 


- 1/2 


+ 1 - 1 - 


7 T~ 


1/4 , -1/4 

m F 0 + m F 0 


(2.46) 
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When tof 0 = 1, these expressions give, 


Tipw 0 (mw 0 = 1) = Trp| 0 (m Fo = 1) = • (2.47) 

in accord with the values predicted in [ 1 ] from the spectral theory/mirror symmetry correspon¬ 
dence. 

We can now use (2.42) to obtain the values of the same traces for local F 2 . We find, 


Trp F2 = 


Trp F2 = 


1 cosh (mp 2 / 2 ) 
4-7T y/mp 2 — 2 

1 


t2 167T 2 
We obtain, in particular 


o cosh 1 (m Fa / 2 ) | 1 \ _ 1 _vr] 

m 2 2 - 4 / m F2 + 2 


1 


TrpF 2 (% 2 = 0) = 


Trp | 2 (m Fa = 0) = ( 3 - 1 ) , 


8 ^ 2 ’ 

•L 

64 \ 7T 


which were already obtained in [14], and 


TrpF 2 (m F2 = !) = —> 


Tr-PF 2 ( m F 2 = 1) = 


1 / 12V3 
432 1 7 r 


- 5 


(2.48) 


(2.49) 


(2.50) 


It can be verified [16] that these values agree with the predictions of the conjecture in [1], 


2.3 Perturbative expansion 

We are now interested in studying the matrix integral (2.31) in the ’t Hooft limit (1.1). As in 
[13], we should first analyze the integrand of (2.31) when h (or equivalently b) is large. At the 
same time, we have to decide what is the appropriate scaling of the parameter appearing in 
the operator, as h becomes large. As it was explained in [13] , in order to recover the topological 
string for arbitrary mass parameter, we have to scale 

logmF 0 ~Ii, h—> 00 . (2.51) 

We recall the variable £ is defined as 


£ 


7 r 

2 h 


log m Fo . 


(2.52) 


This is the mass variable that will be kept fixed in the 
parameter 


1 

g= n ’ 


’t Hooft limit. 


If we introduce the 
(2.53) 
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we can write the matrix integral (2.31) in the form 


Z(N, K) = 


3 -€A/(27rg ) r d N u 

A7! Jr n (2tt) 


JV„, FI (<■)' 4 sinh ( 


JV 


n 

i= 1 




riij2cosh 


where 


'(£) 


V" (it, g) = -glog 

As in [13], we can now use the self-duality of Faddeev’s quantum dilogarithm, 

<h b (z) = $i/ b (z), 

as well as (A. 10), to obtain the following asymptotic expansion for small g, 




Im 


fc >0 


If we write this expansion as 


V{u, g) = Y J Z n vW{i 


£>0 


we find that the leading contribution as g —> 0 is given by the “classical” potential, 
V {0 \u) = — ^ ^ (lmLi 2 (-ie“ +? ) + ImLi 2 (-ie u_? )^ . 


(2.54) 


(2.55) 


(2.56) 


Li 2 _ 2fc (-ie“+«) + Li 2 _ 2fc (-ie“-«)l . (2.57) 


(2.58) 


(2.59) 


The matrix integral (2.54) is an 0(2) matrix model [41], in which the inverse Planck constant 
g plays the role of the string coupling constant, and the potential itself depends on g. In order 
to obtain the ’t Hooft expansion of the free energy, we can use the asymptotic expansion of the 
potential (2.58). In particular, since this expansion only involves even powers of h, we conclude 
that the matrix integral (2.54) admits a standard’t Hooft expansion, of the form 


F(N, h) = log Z(N, h) = Y2 ft 2_2 ^(A, 0, (2-60) 

g> o 

where A is the ’t Hooft parameter introduced in (1.1), and £ was introduced in (2.52). Note 
that, in the planar limit, only the classical part of the potential (2.59) contributes. By using the 
asymptotics of the dilogarithm, one finds that the classical potential behaves as 

V^°\u) « —-, |w| —y oo, (2.61) 

27T 

i.e. it is a linearly confining potential at infinity, similar to the potentials appearing in matrix 
models for Chern-Simons-matter theories [6, 19] and in other matrix integrals associated to 
quantized mirror curves [13]. The potential (2.59), for two values of £, as well as its asymptotic 
form (2.61), are shown in Fig. 1. 

We would like to compute the genus g free energies J r g ( A,£) appearing in the expansion 
(2.60). We will first obtain approximate expressions for the very first free energies, as expansions 
around A = 0, by doing perturbation theory in g, as in [13]. To do this, we regard (2.54) 
as a Gaussian Hermitian matrix model, perturbed by single and double trace operators. The 
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Figure 1. The classical potential (2.59) as a function of u, for £ = 1 (lower line) and £ = 2 (upper line), 
together with their asymptotic form (2.61) when u is large. 


computation is straightforward (see for example [25] for a similar example). For the genus zero 
free energy, we find the following structure, 


W,0 = y (log 


7r 2 Acosh^\ 3 
4 ) ~ 2 


7r 


rim (Li 2 (ie^)) A + ^ / 0)fc A ?l 


k> 3 


-7t(A,£) = -T^log^- -^logA + C'(-l) + 


^(A,0 = 


12 ° 12 
B-2 


k> 1 


2<?(2ff - 2) 


\ 2 - 29 + ^f g , k \ k , g> 2. 


fc>i 


(2.62) 


In writing the second term in the first line, we used the dilogarithm identity 


Li 2 (z) + Li 2 




7) log 2 ( z'j . 


(2.63) 


In the last line, Bo g are Bernoulli numbers. The coefficients f y ^ are themselves non-trivial 
functions of the parameter £. For g = 0, one finds, at the very first orders, 


/o,3 = 7T 
/o,4 = 7T 
/o,5 = 7 t' 


2 1 — 3 cosh(2£) 

24cosh(£) 

4 -73 + 68 cosh(2£) + 45 cosh(4£) 

2304 cosh 2 (£) 

6 534 — 203cosh(2£) — 390cosh(4£) — 165cosh(6£) 
30720 cosh 3 4 (£) 


(2.64) 
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while for g = 1,2, one finds, 


/l,l = ?r : 
/l,2 = 7T' 
/l,3 = 7T 1 
/2,1 = vr' 


2 — 1 + 3 cosh(2£) 

48cosh(£) 

4 127 + 4cosh(2£) — 27cosh(4£) 

2304 cosh 2 (£) 

6 -750 - 265 cosh(2^) + 30 cosh(4£) + 57cosh(6£) 
18432cosh 3 (£) 

6 894 + 577 cosh(2£) + 210 cosh(4£) + 15 cosh(6£) 


(2.65) 


61440 cosh 3 (£) 

These results will be crucial in order to compare the asymptotic evaluation of the fermionic 
spectral traces, to the predictions of [1], 

2.4 The exact planar solution 

The 0(2) matrix model can be solved exactly in the planar limit [17, 42], However, it was noted 
in [19] that instead of using the specific results for the 0(2) case, it is more convenient to first 
consider the 0(n) model for arbitrary n, solve it with the powerful techniques of [18], and then 
take the limit n —> 2. 

In order to proceed, we change variables z = e u in the matrix integral (2.54), and we obtain 


Z(N, K) = 


-++ A 
0 2-rrg^ 


diV - „-i££i(v (o >(*i)+0(g 2 )) IW* - z iY 

F 6 


( 2 . 66 ) 


(2.67) 


m J ( 2 ^) N n ijfa+zj) 

where the classical potential (2.59), when written in terms of z, reads 

y(o) / ) _ log(-g) + ImLi 2 (ize£) + ImLi 2 (ize~£) 

27T 7r 2 

To obtain the planar limit it is enough to consider the classical potential in (2.66). We assume 
that we can model the distribution of eigenvalues by a continuous function on a single connected 
compact support, i.e. we assume that we have a one-cut solution. This is a natural assumption, 
since the potential has a unique minimum at u = 0 and it has a linearly confining behavior 
(2.61). We will take the cut along the segment [a,b] 6 M+. Following the techniques of [18] (in 
the conventions of [19]) we introduce the auxiliary G- functions, 


G<"+) = -i 

a^-\z) = -[ e ^ s+ (z)Gf{z) + e- l ^g + (-z)a^\-z)], 


( 2 . 68 ) 

(2.69) 


where 


<++> 


—1 z 


umi 




sjz 2 - a 2 Vz 2 - b 2 ^ (t 


e ( 2 K 3 with z = asn(u), (2.70) 


9+(z) = 


\Jz 2 — a 2 Vz 2 — b 2 + f-\/e 2 — a 2 Ve 2 — b 2 


z 2 — e 2 


with 


e = asn(i(l — i/)K'). (2.71) 
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Here, K and K 7 are elliptic integrals of the first kind, and , d l (u) are Jacobi theta functions. 
We follow the conventions of [43] for all the elliptic functions and integrals appearing in these 
formulae. The elliptic modulus k and the nome q (= e 17rT ) are given by: 

h = j ) , q = e _7r ^\ (2.72) 

Also, the v parameter is related to the n of the 0(n) model by 

n = 2 cos(7ri/), (2.73) 


so that n —> 2 corresponds to v —> 0. 

Let us denote by C the closed contour encircling the cut at [a, b] clockwise. The following 
equations, due to [18], allow to find the’t Hooft parameter A as a function of the end-points of 
the cut a, b: 


0 = 


1 


2 cos 


dz dW°) 


k(i-v) ) J c 27ri dz 


G {1 ~ u) (z), 


A = 


1 


d -hlWG<% } 


(2.74) 


(2.75) 


2(1 — cos(z^7r)) cos(^) J c 2-k\ dz 

The hrst equation is satisfied if we set b = 1/a, as expected from the symmetry m -H- — rq of the 
matrix integral. So our elliptic modulus is given by k = a 2 . The second equation leads, in the 
limit v —> 0, to the equation 


A = 




7T 


(2.76) 


To determine the function /(£, a), we note that the derivative d/(£,a)/<9£ can be computed by 
deforming the contour and using the residue theorem. Indeed, we have 


d ( dl/(°)' 


dt 


dz 


1 


1 


1 


2-7r 2 i Vz + ie^ z — ie^ z + ie ? 


+ 


1 


— ie~?/ 


z — ie 


(2.77) 


After some calculations, one obtains 
= K 

^ 27ry (a 2 + l) 2 + 4a 2 sinh 2 (£) 


— I Z arcsin 




e^\/T+~a 2 e 2 5 

V a 2 + e 2 £ J \ /a 2 + e 2 ^ 


+ ( Z arcsin 




V a 2 + e _2 £ 


^Vl + a 2 e _2 £ 

Va 2 + e -2 ? 


+ 2a 2 sinh(2£) 


(2.78) 

where Z is the Jacobi Zeta function. The argument of the elliptic functions appearing in this and 
subsequent expressions is now given by the complementary modulus 


k\ = \J 1 — 


a’. 


Since /(£, a) —> 0 when £ —> —oo, we can write 


/(£>«) = J d 


(2.79) 


(2.80) 
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A convenient expression for this integral is in terms of Jacobi theta functions with norne 


One finds, 


/(£,<*)= , lim 

A —>-00 


'1 

rm w ( A ) / 

4 



qi = e < = e mT1 . 


ft 


(2.81) 




where 


+ ^ ^y^a 2 + l) 2 + 4a 2 sinh 2 (£) — \J (a 2 + l) 2 + 4a 2 sinh 2 (A) 
^(0 = F 


(2.82) 


arcsm 


._ , (2.83) 

Va 2 + e 2 «7 

and F is the incomplete elliptic integral of the first kind with modulus k \. This equation deter¬ 
mines the endpoints of the cut as functions of the’t Hooft parameter A. The planar free energy 
is then determined by the equation [18, 19], 


<*F„ , K' 

dV = ~ 2w K =21 ° g?1 ’ 


(2.84) 


up to two integration constants, which can be easily fixed by the weak coupling analysis of the 
previous section. 

The exact planar solution makes it possible to explore the dependence of J-q on the full 
moduli space of A,£. First of all, we can reproduce the perturbative results in (2.62), (2.64) by 
doing a small ’t Hooft coupling expansion. When the ’t Hooft parameter goes to zero, the cut 
collapses to the minimum of the potential. In the z-plane, the endpoint of the cut a goes towards 
1, so we can expand in small k\ = \/l — a 4 . In this case, we can use the qi-expansions of the 
theta functions in (2.82), and after some calculations we find, 


A = 


+ 


1 


-kf + 


1 


, 6 , 115 + 119 cosh(2£) ;8 

' A 9 Q, v 


647T 2 cosh(^)" /i ' 647T 2 cosh(^)' , ' i ' 163847T 2 cosh 3 (£) 

51 + 55cosh(2{) t!0 + o(( , !2) 


(2.85) 


81927T 2 cosh 3 (£) 

Inverting this series and plugging it in (2.84), we obtain, after integrating twice, 


W,£) = co(£) + ( Cl (0 - A + y (log 


7r 2 A cosh£ 


"2 l + 


vr 2 (l - 3cosh(2£)) A s 


24 cosh(£) 


+ 


7r 4 (—73 + 68cosh(2£) + 45cosh(4£)) 4 


( 2 . 86 ) 


2304 cosh 2 (£) 


A 4 + 0(A & ), 


where co,i(£) are integration constants, and we added after integration the missing — £\/(2ir) 
from the prefactor of (2.66). This agrees with the perturbative expansion at genus zero from 
(2.62), (2.64). 

One advantage of the exact solution is that we can also analyze the regime of strong’t Hooft 
coupling. For this, we do an 5-transformation in (2.82) and express our formulae in terms of 


q = e inT = e -i7r / Tl , 


(2.87) 
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which is the relevant variable for the large A expansion. We also do a shift in the integration 
variable to obtain: 


/(£,a) = 



Tp ( 1 — u >( A )/ K ) 
^(l-w( 0 /K) 


dy 


. A /d'i , . i?4 . . 4i?/ N 

U ( » )+ * w+ ^, 


+ ^ ^y^a 2 + l) 2 + 4a 2 sinh 2 (£) — \J(a 2 + l) 2 + 4a 2 sinh 2 (A) ^ 1, (2.88) 


where the elliptic integrals are still evaluated at k\. As we did for the weak coupling expansion, 
we expand the integrand in small q and integrate. After some calculations, we obtain 


1 


k 


1 


^ = log 2 ~7 TX-+ -4 cosh (2^) I 1 - log 7 


1 


87 r 3 


k\ (k 


127T 27T 3 7T 


+ i{ 4 ( 1_l0g l) + 3cosh (46 (-l + 21og^ 



(2.89) 


where we remind that k = a 2 . This can be inverted to yield the series, 


^ _ - 271 V 2 X 
4 


+ 4 + 


rrV A, 


cosh(2£)e 


-47r\/2l 


+ 0 e 




(2.90) 


where we use the shorthand notation 

1 £ 2 
A = vrA + — + 

12 2vr 2 

By using again (2.84), we finally obtain, 

mko = + + ( 2l ® - 1 ) A - 


(2.91) 


1 


4 \ 


4 \ 


32vr 5 


8tt H - — + 7T H - -j= cosh(4£) f e 

V 2 A/ V v 2 A / 


3 -47r\/2l _|_ | e -67rV^A 


(2.92) 


where co,i(£) are integration constants that can be fixed in principle from the weak coupling 
behavior. As anticipated in [1], the strong coupling expansion of the free energy displays the 3/2 
scaling typical of theories of M2 branes [44] , and the coefficient of the leading term agrees with 
the general formula for local del Pezzo Calabi-Yau’s found in [1], The expansion (2.92) is very 
similar to the expansion of the planar free energy of ABJ(M) theory presented in [24]. As we 
will see in the next section, one can in fact recover the result for ABJ(M) theory from (2.92). 
Let us also note that, in the case £ = 0, the formulae above simplify considerably, and one can 
write the periods in terms of indefinite integrals of theta functions, 


A = 27 ri 

= - j dn r 1 $ 2 (2ti) 4 $ 3 (2ti) 2 + c , 


(2.93) 
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where c is again an integration constant. These integrals can be performed in order to obtain an 
expression which is useful for strong coupling expansions, namely, 


A = 


,G 


3,2 


4k 


4vr 4 3 ’ 3 \(k + l) 2 

_ 4k 
dX 167T (k + l) 2 


I I i 

2 ’ 2 ’ 

0 , 0,0 


+ C U 


3 3 4k 

4^3 ( 1,1, —, —; 2,2.2i —- ttt 

4 d 2 2 (A: +1) 2 




4fc 


(A: + l ) 5 


(2.94) 


+ C 2 . 


Another ingredient of the planar solution which can be computed exactly is the density of 
eigenvalues. Let us first consider the resolvent of the 0(2) matrix model (2.66), defined as 


u(p) = — ( Tr—— , , 

N \ p-M / 


(2.95) 


where M is the matrix with eigenvalues Zi, i = 1, • • • ,N, and the bracket denotes the normalized 
vev. We can split this function into its even and odd parts with respect to p, 


u(p) = u+(p) +u-(p). 

The planar limit of the even part can be computed by using the formula [17] 


(2.96) 


(2.97) 


where V ( z ) is the planar potential (2.67), and C is a contour around the cut [a, 1/a] anti-clockwise. 
We find, 


<4 (p) = -TW 


i r dw v'(w)w \/(p 2 - a2 )(p 2 - 


2A Jc 2 tt\ p 2 - w 2 ! {W 2 _ a 2 ){w 2 _ A_ 2) ' 


w+(p) = - 


a-\J(p 2 ~ a 2 ){p 2 - 4y) 


47r 2 ip 2 A 


n 


a . ie s 
n , arcsm — 
p z a 


a 4 +n 


a . le 
, arcsm — 
p z a 


-e 


° 4 - L(p , a) 


+ 


1 


4tt 2 pX 

where II is the elliptic integral of the third kind, and 


( arctan-arctane^p ) , (2.98) 

VP J 


L(p,a) = lim II ( —,iA 




P 


a 4 I = 


-l 


a 2 (a 2 — p 2 ) 


nll -!V 


i - 4 - k' 4 


When z E [a, 1/a], the eigenvalue density is given by the discontinuity equation 
p(z) = 4 (z + iO) - u+tz - iO)) 


(2.99) 


( 2 . 100 ) 


ay / (z 2 - a 2 )(4 - z 2 ) 


27r 3 iz 2 A 


a 2 ie£ 

II I —. arcsin — 


a 4 ) + II —, arcsin 


le 




a 4 j — L(z, a) > . 


( 2 . 101 ) 


This expression can be checked by doing a numerical simulation of 300 eigenvalues relaxed 
into a configuration which approximately minimizes the effective action of the matrix model, 
as in [45]. The results are shown in Fig. 2, after going back to the initial u variable, so that 
p(u) = p(z(u))dz/du, with z(u) = e u . 
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Figure 2. Left: eigenvalue density p(u) for A = 1, £ = 0 against a histogram showing the numerical 
density of IV = 300 relaxed eigenvalues. Right: same plot with A = 3/4 and £ = 5. 


2.5 Relation to the ABJ(M) matrix model 

In [46], a matrix model computing the partition function of ABJ(M) theory [47, 48] on § 3 was 
derived, by using localization. This model turns out to be closely related to the topological string 
on local P 1 xP 1 , In the case of ABJM theory, this was first noted and exploited in [6, 24] in the 
context of the ’t Hooft expansion, and then in [9, 10] for the M-theory expansion (which involves 
as well the refined topological string in the NS limit). The generalization to ABJ theory was done 
in [23, 49]. Since the matrix model (2.31) gives a non-perturbative completion of the partition 
function on this geometry, it is natural to wonder whether it is related to the ABJ(M) matrix 
model. In fact, one can recover the ABJ(M) matrix model from (2.31) provided one considers 
complex values of the parameter m^ 0 . Note that the operator (2.10) is no longer self-adjoint in 
this case, and in addition one has to be careful with the resulting multi-valued structure, since 
the integral kernel depends on the logarithm of m^ 0 . 

Let us then set 


log mw 0 = i h — 2vri M, M e Z>o- (2.102) 

Here, the integer M will be identified with the difference between the ranks of the two Chern- 
Simons theories in ABJ theory [48] . This relationship is the one suggested by the explicit results 
of [23, 49]. In these papers, the grand potential of ABJ(M) theory is written in terms of the 
topological string on local Po- If we now use the explicit expression (2.20) for the integral kernel, 
we find 


PW o 


_ g—i7t/4+i7rM/2g\/irftq/2 



+ /h) i 

— 2 cosh(\/7r h p) 


$ v /n/L(9 + ^My/Trjh- \\fhjn) 


(2.103) 


Due to the form of the arguments, we can use the functional equations for the quantum diloga- 
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rithm, (A.9a) and (A.9b), and we obtain 


PFq = e -ifc/4+i7rM/2 e V^q/2 


/ M—l 

2 1 

TT 1 

__M+1 1 _|_ ^(q^/Vfh+ms/K) I 2cOsh(\/7T^p) 

V s - 


M — l 
2 


(2.104) 


\ s = 


_|_ e 2n(qy/7r7H+i7rs/h+i/2)j _ 

_ A/f_LI I 1 


_ -M+l 
2 


_|_ ]_jMg2\/7rSq 


^Vwhq/2 


To make contact with ABJ(M) theory, let us define 

h 


k = 


IT 


(2.105) 


and let us introduce the variables 


u = 27r\/fcq, y = 2TrVk\ 


(2.106) 


so that [u,v] = 2TT\k. In these new variables, and after a similarity transformation, we find, 


where 


Pabj(m) 


Ap ¥o A -1 =e -™k/i+mM/2 


2cosh(v/2) e t + (_i)M e -i 


2 

Pabj(m)? 

M — l 


(2.107) 

|j tanh ( 

-M+1 

< u + 27ris^ 

(2.108) 

v 2 k ) 


is, up to a similarity transformation, the operator appearing in the Fermi gas formulation of 
ABJM theory [6] and of ABJ theory [21-23]. Since the phase appearing in (2.107) is the same 
one appearing in the relation between Fq and F 2 , we also conclude that, 


Pf 2 — Pabj(M)> (2.109) 

up to a combination of unitary and similarity transformations. The dictionary between the 
parameters is (2.105) and 

/irk \ 

7T7.f 2 = 2 cos ( —- irM J . (2.110) 

In particular, the spectral traces of the kernel of the ABJ(M) matrix model can be obtained 
from the traces of the F 2 operator. This can be easily tested for ABJM theory, in the case 
k = 2, M = 0, by using the expressions (2.48). The relevant value of the mass parameter is 
mp 2 = —2, which is a branch point for the functions in (2.48). However, the traces at this point 
are well-defined and one finds the correct values [7] 

Tr pabjm | fc=2 = g. TtPabjmL^ = g4 “ ^2' ( 2 ' m ) 

We also find that, when k = 2 and M = 1, which is the maximally supersymmetric ABJ theory, 
the theory is equivalent (at the level of spectral traces) to the maximally supersymmetric case 
of local F 2 with m ¥2 = 2, or equivalently of local Fq with m ¥o = 1. 
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The exact results for the planar solution found in the previous section can be used to re¬ 
derive the exact planar solution of the ABJ(M) matrix model, first obtained in [24, 50]. Indeed, 
due to (2.102), the exact planar free energy of ABJ(M) theory can be obtained from the general 
formulae obtained above by setting 


i7r m 2 M 

~2 n~ 


( 2 . 112 ) 


As a check, note that the shifted variable (2.91) becomes 




i 

24’ 


(2.113) 


where 



M 

k 


(2.114) 


has to be identified as the B field of ABJ theory. The shift (2.113) is precisely the one found in 
[24], In addition, the strong coupling expansion (2.92) becomes, 


M\P) = 


\/2 A 3 / 2 

37r 


+ Cq + 




P + P~\-2,V7X 

8t r 4 


- ^4 | ^(/ 3 2 + 16 + ^" 2) + + ^ _1 ) 2 } e_47a/21 + ° , 

where 

/3 = e ~ 2 * iM / k . (2.115) 

The function in (2.115) is precisely —l/(47r 4 ) times the planar free energy obtained in [24]. 

This overall factor is due to our different conventions for the string coupling constant. 

It should be noted however that the perturbative and weak coupling expansion worked out 
for the matrix model (2.54) can not be used for ABJM theory in the form presented above. 
For ABJM theory, M = 0, so that £ = i7r/2, and the expansions (2.64), (2.65) and (2.86) 
diverge. This is not a problem of our exact solution, but rather a breakdown of the Gaussian 
approximation. The reason is that, when considering the particular limit of ABJM theory, the 
“classical” potential (2.59) is no longer a perturbed Gaussian, since it is exactly given by the 
r.h.s. of (2.61), and in particular it is not smooth at u = 0. One can however still obtain the 
correct weak coupling expansion from the exact planar solution, and one obtains, 


A = 


8^ 2 mi+ 16^2 


m\ + 


65 


1536tt 2 


m\ + £>(m 4 ), 


where mi is given by (2.79), as well as 


(2.116) 


-To (A) — cq + ci A + 



3\ x2 7T 4 x4 28377% 6 96l7T 12 'g 
2 ) 9 A + 5400 A 19845 A 


which is precisely (up to an overall factor — 1/(47 t 4 )) the expression found in 
one finds the relations 

A = -^-r J dri 7? 3 (2ri) 4 i? 2 (2Ti) 2 , 

= - J dn 7~i i? 3 (2ti) 4 i? 2 (2ti) 2 + c, 


C>(A 10 ), (2.117) 

[24], In addition, 


(2.118) 


which are obtained from (2.93) by exchanging i? 2 with This can be also integrated explicitly, 
as in (2.94), and the result is in precise agreement with the result of [24], 
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3 Comparing the matrix model to the topological string 

3.1 Predictions from the spectral theory/mirror symmetry correspondence 

The conjecture of [1] gives a very precise prediction for the ’t Hooft expansion (1.2) of the 
fermionic traces of the trace class operators obtained by quantizing mirror curves. We will now 
summarize some of the results of [1], specialized to the case of interest, namely local P 1 x P 1 (see 
also [13, 14] for other summaries of the main results of [1]). According to the conjecture of [1], 
the basic quantity determining the spectral properties of the operator pf 0 is the modified grand 
potential J(/x,m ¥o ,h). This function depends on the “chemical potential” /i, which is related to 
the “fugacity” k entering in (2.26) as 

K = e (3.1) 

as well as on the mass parameter mp 0 appearing in (2.10). The modified grand potential is 
determined by the enumerative geometry of the CY. We first need a dictionary between the 
parameters ix, m Fo , and the parameters appearing in the enumerative geometry of local P 1 x P 1 . 
This CY has a “diagonal” Kahler parameter T, which is related to /j by 


T = 2/igff. 


(3.2) 


Here, the “effective” /i parameter is determined by the so-called quantum mirror map of [4] (see 
[ 1 ] for the notation), 

/Jeff = fx~ ^ ^a^(^)e -2 ^. (3.3) 

Z t>4 

In this paper we will not need the explicit expression of this quantum mirror map, since it is 
not relevant in the ’t Hooft limit we will focus on. In addition, there is a Kahler parameter 
T m associated to the mass parameter Geometrically it measures, roughly speaking, the 

difference in sizes between the two spheres in local P 1 x P 1 , and one has the relation 

T m = — log m ¥o . (3.4) 

The modified grand potential has the form, 

J(lx,m ¥o ,h) = j( p \fx eS ,m ¥o ,h) + JM 2 (Heft,m Fo ,h) + Jws(/4ff, rn ¥o , K). (3.5) 

Here, J^ p \fx, m ¥o , h) is the perturbative part, which is a cubic polynomial in /j: 

_ 1/YV+ (Y - Y) it + Afa.h). (3.6) 


When m ¥o = 1 one recovers the expression presented in [1]. When m ¥o 7 ^ 1, the part which 
depends on [x can be obtained in a relatively straightforward way by working out the classical 
grand potential [16, 51], or by using the result for local P 2 and the dictionary between this model 
and local Po- A precise expression for the function A(m ¥o ,h) has been obtained by Y. Hatsuda 
[51]. His expression exploits the relationship between ABJ theory and topological string theory 
on local P 1 X P 1 discussed in section 2.5. It is given by, 


A(m Vo ,h) 


log 3 m Wo log m Vo / 7 t _/i_\ fh 

A&'kH 4 V3/i 127T/ c l 7r 



(3.7) 
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Let us spell out the details of this formula. The function A c (k ) was introduced in [6] in the Fermi 
gas approach to ABJM theory. It was determined explicitly in [52] and further simplified in [53]. 
It reads, 

= + P-8) 

In (3.7), Fcs(k,M ) is an analytic continuation of the Chern-Simons free energy on the three- 
sphere for gauge group U(M) and level k, 

Fcs(k, M) = log Zcs(k, M ), (3.9) 


where M is related to the parameters of our problem as 


M = 


h + i log m Fo 
2tt 


(3.10) 


Note that this is precisely the relation (2.102) used in section 2.5. As it is well-known, the 
Chern-Simons partition function for integer M is given by [54] 


M 


Zcs(k, M) = & _M / 2 JJ ( 2 sin 
3 = 1 


7TJ 

k 


M-j 


(3.11) 


but in view of (3.10) we have to extend it to arbitrary complex M. This can be done in various 
equivalent ways [51, 55, 56], but in this paper we will not need the precise form of this extension. 

The “membrane” part of the potential JM 2 (/Aff> w-F 0 , K) appearing in (3.5) will not be relevant 
for our purposes. It is fully determined by the refined BPS invariants of the topological string 
in this CY background, see [1, 13] for details. Finally, the worldsheet part of the modified grand 
potential is 

Jwsik'es, rn ¥o ,h) = e~^ vd ' T , (3.12) 

g>0u=l d V ' ' 

where 

T = (T,T m ), (3.13) 

n d are the Gopakumar-Vafa invariants [57] of local P 1 xP 1 , and d = (di, efe) are tw° non-negative 
integers. 

One of the consequences of the conjecture of [1] is that the fermionic spectral traces Zp 0 (IV, h) 
can be obtained as integral transforms of the modified grand potential, 

Z WQ (N,h) = — [ (3.14) 

2vri J c 

where the contour C goes from e -17r/,3 oo to e 17r / 3 oo (so that the integral is absolutely convergent). 
The formula (3.14) leads to a precise prediction for the ’t Hooft limit of the fermionic spectral 
traces. Note that, in order to keep the dependence on both Kahler parameters, we have to take 
a limit in which logmif 0 scales with h , as we required in (2.51). We then consider the ’t Hooft 
limit of in which 


and 


fi oo, 

= C fixed, 
h 


771]Po —^ OO, 

h —> oo, 

(3.15) 

i logmF ° 

= £ fixed. 

(3.16) 
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The parameter £ was introduced in (2.52). We will express often the results in terms of 


m = e 


4£ 


(3.17) 


which corresponds to the mass parameter appearing in the standard topological string free en¬ 
ergies. Indeed, with this definition, one has that 

27T 

- log(m) = —T m . (3.18) 

n 

In the ’t Hooft limit, the membrane part of the grand potential in (3.5) goes to zero, and 
/r e ff —> g- The remaining ingredients have non-trivial’t Hooft-like expansions. The expansion of 
A(rri ¥ Q , h) can be easily worked out. The function A c (k) has the large k expansion [52]: 

M k ) = -g^C(3) + 2 log(2) + 2C'(-1) + 6 (^) + E W ^(-lr 1 ^, (3.19) 

^ <7>2 V J 


where 

BonBon—2 

° 9= (Ag)(2g-2)(2g-2)V 

On the other hand, in the limit we are considering, M —> oo but 


27T 2 i 

h 


M = 7ri — 


2£ 


(3.20) 


(3.21) 


is fixed. This is the standard’t Hooft expansion of Tcs(^/ 7r ) M), worked out at all genus in [26], 
and with’t Hooft parameter (3.21). One then finds an expansion of the form, 


A(m Wo ,h) = ^2A g (£)h 2 2ff , 
g> o 

where Hi(£) includes as well a logarithmic dependence on h , and 
C(3) - 2Li 3 (-e 2 «) 


ako = 


87T 4 


A(£) = ^ log (I6vr 2 cosh^) - ^ log H + C'(-l), 


A g (0 = (2vt 2 ) 2 5- 2 (-1)s- 1 |(49 - 2)c 9 - 


B 


2-7 


2g{2g-2)\ 


Ll3-2g ( —e 


-P 2 € 


It follows from this expression that 

4,(0) = 2(-l y- 1 (4tt 2 ) 29 “ 2 c g (3 - 2 3 - 2 ^) , g > 2, 
in agreement with the result presented in [1] for m Fo = 1, 


A (m Fo = 1 ,h) = -A c - - A — 


h 


7 r 


2 h 


TV 


(3.22) 


(3.23) 


(3.24) 


(3.25) 


One concludes that, in the’t Hooft limit (3.15), the modified grand potential has the asymp¬ 
totic expansion, 


j’ tHooft (c,£y = Ea(caa 2 “ 29 , 

<7=0 


(3.26) 
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where 


J »K-« = l < 3 - ~ ++ 

MCO = |c + ^i(e) + ^ nst (i,m), (3.27) 

jg(CO = MO + (47 t 2 ) 29 ~ 2 F^ (t,m ), g> 2. 

Here, we have introduced the variable 

t = 4ttC, (3.28) 

and F™ st (t, rri) is the worldsheet instanton part of the standard genus g topological string free 

energy. In order to obtain the ’t Hooft expansion of the fermionic trace and make contact with 
(2.60), we have to calculate the integral in (3.14), by doing a saddle-point expansion for h large. 
Let us denote by 

^J-J°(A,m)^ 2 -^ (3.29) 

9> 0 

the asymptotic expansion of the logarithm of the integral in (3.14). At leading order, one finds 


_ dJ 0 ((,m) 
d( 


(3.30) 


which determines the ’t Hooft parameter A as a function of £, and conversely, ( as a function of 
A. The genus zero free energy J 7 ^ 0 (A, rn) is then given by a Legendre transform, 

^o°(A ,m) = J 0 (C(A),m) - AC(A). (3.31) 


In particular 


_ 

d\ 


The next-to-leading order correction to the saddle-point, .Ff 0 (A, m), is given by, 

M (A, rn) = Jr(C(A),m) - ^log ^2vr^l^ . 


(3.32) 


(3.33) 


The higher order corrections can be computed systematically by using the results of [58] (already 
exploited in this context in [13]): in the saddle point approximation, the integral in (3.14) 
implements a transformation from the large radius frame, appropriate to J g ((, m), to the conifold 
frame. Therefore, the functions J r J°(A,m) appearing in (3.29) are the topological string free 
energies of local P 1 x P 1 in the conifold frame. According to the conjecture of [1], they should 
be equal to the matrix model free energies appearing in (2.60) 

F g (\,m) = jJ 0 (A,m), g > 0. (3.34) 

This was tested in detail for local P 2 and for local P 2 (for a fixed valued of mf 2 ) in [13]. We will 
devote the rest of this paper to an explicit verification of (3.34), and in the next section we will 
compute the r.h.s. of (3.34) by using standard techniques in topological string theory. 
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3.2 Topological strings on local P 1 xP 1 

Let us review some basic facts about the special geometry of local P 1 xP 1 . Since this has two two- 
cycles, we can regard it as a two parameter model, and its mirror has then two complex moduli 
z\, Z 2 ■ However, it has been known for some time that m = z\lz 2 does not receive quantum 
corrections, therefore it should be rather regarded as a parameter (see for example [30]). We will 
then have a complex modulus, z = z 2 , and a “mass” parameter m. The periods will be obtained 
as solutions to a single Picard-Fuchs (PF) equation corresponding to the operator [30]: 

C = ^8(1 — m) 2 z 2 — 4(1 + m)z + + ^16(1 — m) 2 z 2 — 4(1 + m)z^6 2 

+ ^6(1 — m) 2 z 2 — (1 + m)z\6, (3.35) 

where 

9 = z —. (3.36) 

az 

This is the form of the operator which is appropriate for the large radius point z = 0. As usual 
in local mirror symmetry, there will be a constant solution 1, a logarithmic solution 

gi(z) = log(z) + 01 (z), (3.37) 


and a double logarithmic solution, 

g 2 (z) = log 2 (z) + 21og(z)oi(z) + 0 2 (z). (3.38) 

In these equations, 01 , 2 ( 2 ) are power series around z = 0, whose coefficients depend on m. The 
very first orders read, 

01 (z) = 2(7772 4 - m~2)m2z + 3^(m + mF l ) + Ajmz 2 

20/. 3 _ 3 . 1 _ 1 . \ 3q . A 

H-(m 2 + m 2 ) 4 - 9(7772 + m z + 0(z ), 

3 i (3.39) 

02(2) = 4(m2 4- 7 ? 7 _ 2 )??7 2 z + ^13(777 4- mnT 1 ) + 40^?77Z 2 

+ ^ ^41(?772 4- 777 , — 2 ) 4- 279(7772 + 777 _ 2 )^ 777 2 Z 3 4- 0(z 4 ) 

Let us now consider the following linear combinations of the basic periods, 



(3.40) 


(3.41) 


The first A-period determines the flat coordinate t through the mirror map, while the second, 
B-period determines the genus zero free energy Fo(t,m) at large radius, 


-t = n 


(lr) 
A > 


dF 0 

dt 


— n (lr) 

— ri B . 


(3.42) 
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After integration, we get, 


t 6 log m o i _i i _ t 1 / _i \ 

Fo(t,m) = —-—— t — 2 (ms +ra 2 jm 2 e —+ ) + 16Jme 

/ 3 g l\3 . 

: ^(ra! + m _ 2 ) + 81 (m 2 + m _ 2 )Jm 2 e~ 3t + 0 (e~ 4t ) 


2 

27' 


Equivalently, one can obtain the same information from the equations 

dt 2 ^ ^ lby^z ^ 

<9^ ttz^I -4( v /m + l) 2 z V4(Vm +l) 2 z - lj ’ 

<9 2 F 0 __2_ K / 4(y / m + l) 2 z - 1 \ 

dtdz Zy /l - 4(y/m - l) 2 z V 4 (\/^— l) 2 ^r — !/ 


(3.43) 

(3.44) 


(3.45) 

(3.46) 


which can be obtained from the results of [59] for local F 2 , together with the dictionary relating 
the moduli of local F 2 to those of local Fo- 

We now analyze the theory near the conifold locus given by the vanishing of the discriminant: 


A = 1 — 8 (m + 1 )z + 16(m — l) 2 z 2 = ^4(1 + \fm) 2 z — 1 j ^4(1 — \fm) 2 z — 1 j. (3.47) 


Note that there are two different branches of the conifold locus, related to the two square roots of 
m. For each value of m, we have a different conifold point in each of the branches of the conifold 
locus, and we have to analyze the topological string near an arbitrary point, as a function of m. 
For m = 1, the topological string near the corresponding conifold point at z = 1/16 has been 
analyzed in [13, 24, 60]. We will pick for convenience the branch of positive roots, and introduce 
the local variable, 

y = 1 — 4(1 + ^/m) 2 z, (3.48) 

which vanishes at the conifold point 


1 

Zc_ 4(1 + \/m ) 2 ' 


(3.49) 


In the variables appropriate to the conifold point, the PF operator becomes 
C = 4(y-l) 2 [y(y-l) 2 + 4y)d 3 y 

+ 4(y - 1)(2 y 2 (y - l ) 2 + y( 1 + y) 2 + 8 ^0 2 y (3.50) 

+ ^3 y 3 {y — l ) 2 + 4 y 2 y + y (+ 1) + 16 y^jd y . 


There is a basis of solutions given by a constant solution 1, a vanishing solution 


h(y) = y + 0(y 2 ) 


(3.51) 


and a logarithmic solution 


f 2 (y) = log(y)/i(y) + s(y), 


s(y) = y + 0{y 2 ). 


(3.52) 











It is easy to solve for f\ (y) and s(y) as power series in y with m-dependent coefficients: 


h(y) = y- 
s{y ) = y - 


cosh(2£) — 11 9 9cosh(4£) — 124cosh(2£) + 827 3 4 

16 r+ 1536 V 

7cosh(2£) - 45 2 27cosh(4^) - 380cosh(2£) + 1561 3 4 

32 V + 1152 V + ^ ' 


(3.53) 


where we expressed the results in terms of the variable £, related to m by (3.17). The analytic 
continuation of the large radius periods to the conifold point must be a linear combination of the 
two solutions fi(y), fiiy) found above. By expanding the exact results (3.45)-(3.46) around the 
conifold locus, one finds 


n? } (*) = 




n£ } (z) = 



7T cosh £ 




(3.54) 


(3.55) 


where C±, C 2 are a priori ^-dependent constants which we do not know how to evaluate analyt¬ 
ically. These constants are given by the values of the large radius periods at the conifold point, 
i.e. 

C 1 = n%\z c ), C 2 = U { ^(z c ). (3.56) 

The constant C\ can be computed analytically, and we present this computation in Appendix 
B. The constant C 2 can be calculated numerically, by evaluating the series (3.39) at the conifold 
point (where the series still converges). However, as we will see in the next section, the value of 
these constants is predicted by the conjecture of [1], and we will find a precise agreement with 
the analytical and numerical evaluations of Ci 2 , respectively. 

The above results determine the genus zero free energy at large radius Fo(t,m). The genus 
one free energy can be obtained, for example, from the result for local F 2 in [59], by using the 
map of moduli. One finds, 

Fi(t,m) = log (n2^z 7 (16(m - 1 ) 2 z 2 - 8(m + 1 )z + 1)) - ^ log > (3.57) 

with the large radius expansion 

Fi(t,m) = — log (m) + ^ - ^(1 + m)e~ t + ■■■ (3.58) 

The higher genus free energies near the conifold point can be obtained by integrating the holomor- 
phic anomaly equation. A systematic computation has only been performed at m = 1 [24, 60], 
but it will provide us with useful tests, as we will see in the next section. 


3.3 Comparison 

Let us now compare the results from the matrix model/spectral theory side, with the predictions 
from the conjecture of [1]. We start with the genus zero free energy. By using the expansion 
(3.43), we can write 

J ° = 4^) = 16^4 (^’™) “ + 167r4j V£)l > ( 3 - 59 ) 
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where t is related to C, by (3.28) and it is the standard Kahler parameter of the geometry. We 
also recall that m is related to £ by (3.17). The ’t Hooft parameter is given by (3.30), which 
reads in this case, 

47T 3 A = < r) -y. (3.60) 

The form of the matrix model expansion (2.62) indicates that A must be a vanishing period at 
the conifold point, i.e. 

n ^\z = z c ) = C 2 = ^. (3.61) 

This gives a prediction for the value of the constant C 2 , which we have verified by evaluating 
this constant numerically. This test involves doing a high precision numerical sum of the large 
radius expansion of II^, 1 '*( 2 : = z c ), for different values of m. 

We should note that, as far as we know, the constants that have to be added to the H-periods 
in order to obtain a vanishing period at the conifold point are not known a priori, and they have 
to be determined on a case-by-case basis, often numerically (see [61, 62] for some examples.) 
According to the conjecture of [1], this constant is obtained from the terms in the modified grand 
potential which are linear in the chemical potential and next-to-leading in h. On the other hand, 
it follows from the general construction in [1] that these terms are in turn determined by the 
linear terms in the moduli appearing in the large radius free energy T)^. Therefore, we have the 
following consequence of the conjecture of [1]: in a toric CY threefold, the constant terms which 
appear in the vanishing periods at the conifold are determined by the coefficients of the linear 
terms in F^. This is of course also the case for the example of local P 2 studied in [13], and 
provides an intriguing link between refined genus one free energies and the special geometry of 
the conifold point. 

Let us now proceed to compute from (3.32), which reads 


9^0” = ± n m 

d\ 4 tt a ' 

This can be integrated to find, up to an integration constant, 


(3.62) 




A 2 


7r 2 Acosh£\ 


2 VV 4 ) 2/ + 4tt 

4 —73 + 68cosh(2£) + 45cosh(4£) 4 


C\(0 x , _ 2 1 - 3cosh(2£) A 3 


+ 7T 


A + 7T 


24 cosh(£) 


-A 


+ 7T 


2304 cosh 2 (0 

6 534 — 203cosh(2£) — 390cosh(4£) — 165cosh(6£) 5 


30720 cosh 3 (^ 


A 5 + 0(A b ). 


(3.63) 


The A-independent function Ci(£) appearing here is the one appearing in (3.55). This result 
agrees with the results in (2.62), (2.64), and (2.86) obtained in the matrix model, provided that 


Ci(£) = ”Im (Li 2 (ie 5 )) . (3.64) 

Note that, as in the examples of [13], the r.h.s. of (3.64) is a prediction of the conjecture of [1] 
on the value of the A-period at the conifold point. This prediction comes from the explicit form 
of the potential (2.59), which is in turn determined by the explicit form of the integral kernel 
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(2.22). As shown in Appendix B, the value (3.64) agrees precisely with the analytic evaluation 
of the A-period at the conifold. 

Finally, we note that, according to the result (2.62) in the matrix model, the integration 
constant in J-jJ 0 (A, rri) should vanish. This implies that 


0 = ^°(A = 0,m) 


(yo(C(A),m) - AC(A,m)) 


A=0 


(3.65) 


Using (3.59), this yields the following relation between the function Ao(£) in (3.23) and the value 
of the genus zero free energy at the conifold point: 


A)(0 





(3.66) 


where we used that t(z c ) = — Ci(£). This is yet another remarkable consequence of the conjecture 
of [1] for the special geometry of the conifold point, and we have verified it numerically with high 
precision. To give a flavor of the validity of (3.66), in Fig. 3 we show the value of Ao(£), as given 
in (3.23), against a numerical evaluation of the r.h.s. of (3.66) for some values of £. 



Figure 3. The continuous line shows the exact function Ao(£), as given in (3.23), while the dots are 
numerical evaluations of the r.h.s. of (3.66) for some values of £. 


We conclude that, at the level of the genus zero free energy, there is a remarkable agreement 
between the result obtained from the matrix model (i.e. from the spectral theory side) and the 
predictions of [1] based on topological string theory. In particular, one can use the results of the 
matrix model/spectral theory side to obtain non-trivial information about the conifold theory 
which is not available otherwise (as far as we know): the analytic values of the periods and the 
genus zero free energy at the conifold point. 

Let us now consider the genus one free energy. Note that 


MC,0 = ^C + Ai(0 + Fr\t,0 = Ai(0 + | + (3.67) 

3 6 


where we used (3.57) and 
equation 


the expansion (3.58). If we now 


d 2 MC) 

d( 2 


A dX 
— 47T- 

dt : 


take into account (3.33) and the 

(3.68) 
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we obtain 


^f°(A) = ^i(£) + | - ~ logm - ^log (m 2 z 7 (\){16(m - l)V(A) - 8(m + l)^(A) + 1)) 

1 / 1 dz( A)\ 

+ 2 bS ( “ 8^ d\ ) ' 

(3.69) 

By using the explicit expression for Hi(£) in (3.23), we find that the small A expansion of this 
function is, 


•7?°(a) = -; o io g A- ; logh+c^-i) 


— log A —- 
12 6 12 

2 —1 + 3cosh(2£) 4 127 + 4cosh(2£) — 27cosh(4£) 2 


(3.70) 


+ vr 


48 cosh 2 (£) 


-A + vr" 


2304 cosh 2 (0 


A 2 + 0(A 3 ), 


which is in precise agreement with what was found in (2.62), (2.65). 

There are some non-trivial tests that can be done at higher genus, in the case £ = 0, by 
using the results in [24, 60]. Let us first recall that the topological string free energies in the 
conifold frame, when expanded around the conifold point in terms of a vanishing period, have 
a universal critical behavior characterized by a pole of order 2 g — 2, for g > 2 [63]. It was then 
pointed out in [64] that the full expansion satisfies a “gap condition,” i.e. after this pole, the 
rest of the expansion is regular and it starts at zeroth order in the vanishing period. This has 
been exploited to constrain solutions to the holomorphic anomaly equations. However, in the 
matrix model free energies, as one can see in (2.62), the expansion around the conifold fulfills a 
“strong” gap condition, in the sense that the expansion in A after the pole starts at first order in 
A (and not at zeroth order). In contrast, the conventional topological string free energies satisfy 
only a “weak” gap condition. In practice, this has the following consequence. Let us consider 
the instanton part of the large radius, genus g free energies Fg nst (t,m), and let us perform a 
symplectic transformation to the conifold frame. The “weak” gap condition of [64] implies that 
the expansion of the resulting quantities around the conifold point is of the form, 


B 


2 g 


.2-2g 


2 ^( 2 g - 2 ) 


bg(0 + 0(t c 


(3.71) 


where t c = 47 t 2 A is a vanishing period at the conifold 3 . Then, it follows from the last line in 
(3.27) that 

0 = 2g(2g-2) X2 ~ 29 + (^^MO + MO + ( 3 - 72 ) 

Therefore, consistency with the expansion (2.62), which satisfies a strong gap condition, requires 
that 

= 9 >2. (3.73) 

This can be regarded as yet another prediction of spectral theory for the topological string, since 
the coefficients A g (£) have been fixed by consistency with studies of the spectrum. For £ = 0, 

,5 We are considering just the instanton part of the large radius free energies, so we are not including the constant 
map contribution to these amplitudes to zero. Note however that adding this contribution does not lead in general 
to a strong gap condition at the conifold. In other words, — b g is not the constant map contribution at large radius. 
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the constants b g (0 ) can be computed systematically from the holomorphic anomaly equations 
[24, 60]. One finds, for the very first genera, 


b 2 (0) 


1 

1152’ 


63 ( 0 ) 


23 

5806080’ 


64 ^ 0) 278691840’ 


(3.74) 


and by using (3.24), one verifies that (3.73) is indeed satisfied. 

Finally, we note that the genus two free energy in the conifold frame is given by [24, 60], 


F? st (t c ,a = 0 ) 


1 1 | 53 t c 2221f 2 

240 t 2 c ~ 1152 + 122880 ~ 14745600 + ' " 


(3.75) 


The third term in this expansion agrees with the coefficient / 2.1 in (2.65), for £ = 0 after taking 
into account the overall factor (47 t 2 ) 2 in (3.27). We conclude that the’t Hooft expansion of the 
fermionic traces, as calculated by the matrix model, is in perfect agreement with the predictions 
of [1] (and with the result of [51] for the function A(m® 0 ,K)). 


4 Conclusions and open problems 

In this paper we have found an explicit expression for the integral kernel of the trace class 
operator associated to the mirror curve of local Fo, for arbitrary values of the mass parameter. 
This makes it possible to obtain a matrix model computing the fermionic spectral traces of this 
operator. This model turns out to be an 0(2) model, which can be exactly solved in the planar 
limit. Using this matrix model, we have verified in detail that the fermionic spectral traces of 
(2.20) provide a non-perturbative definition of the topological string on this geometry, in the 
sense that their asymptotic’t Hooft expansion is given by the genus expansion of the topological 
string. This provides yet another test of the conjecture in [1], In particular, our calculation checks 
the conjectural form of the quantum-mechanical instanton corrections to the spectral problem 
proposed in [1]. 

There are various obvious problems raised by our results. It would be interesting to improve 
our checks by calculating higher genus amplitudes directly in the matrix model, although this 
type of calculations are not simple for 0(n ) models. Even at genus zero, it would be interesting 
to have an analytic proof that the function J -0 obtained in the matrix model agrees exactly with 
the genus zero free energy of the topological string J r g°. Another obvious open problem is to 
extend this type of calculations to other geometries, like for example local B n , where B n is the 
blow-up of P 2 at n points. To do this, we would need an explicit form for the integral kernels of 
the corresponding trace class operators. It would be also interesting to find exponentially small 
corrections to the ’t Hooft expansion of the matrix model studied here, in order to construct a 
trans-series expansion of the matrix model free energy, in the spirit of [65] (see [66] for a recent, 
detailed case study of trans-series in the quartic matrix model). This could then be compared 
to the predictions of [1] and/or to the trans-series construction of [67, 68]. 

Another research direction concerns the field theory limit of the model analyzed in this 
paper. It can be explicitly shown that the spectral problem for the operator (2.10) has a double¬ 
scaling limit in which one recovers the quantum spectral curve of 51/(2) Toda given in [69]. This 
corresponds to the held theory limit of the topological string, which is pure N = 2 Yang-Mills 
theory [32]. According to [3], the NS limit of the instanton partition function of [70] provides a 
quantization condition for this spectral problem. This can be verified by using the perturbative 
WKB approach [5], but there are also non-perturbative corrections (see [71-73] for different 
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perspectives on this issue). It would be interesting to analyze this field theory limit by using the 


tools introduced here. 

Finally, as we have explained, the matrix model in this paper generalizes the ABJ(M) matrix 
model, and in particular extends it to arbitrary values of M. This is due to the fact that the 
dependence on M is through the mass parameter mf 0 , as shown in (2.102). In contrast, in the 
existing matrix models for ABJ theory [21, 46], M has to be in principle a positive integer. This 
might be useful in order to relate ABJ theory to higher spin theories [74, 75]. 
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A The quantum dilogarithm 

The quantum dilogarithm 4>b(x) is defined by [33-35] 



(A.l) 


where 




, Im(b 2 ) > 0 


(A.2) 


and 


c b = ^ (b-f b x ) . 


(A.3) 


An integral representation in the strip |Imz| < |Imcb| is given by 



(A.4) 


Remarkably, this function admits an extension to all values of b with b 2 0 R<o- 4 ) b(®) is a 
meromorphic function of x with 


poles: Cb + iNb + iNb 1 , zeros: — Cb — iNb — iNb 1 . 


(A.5) 


The functional equation 



*b(0)= (% 
\Q 



(A.6) 


allows us to move 4>b(x) from the denominator to the numerator. In addition, when b is either 
real or on the unit circle, we have the unitarity relation 



(A.7) 
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The asymptotics of the quantum dilogarithm are given by [76] 


$b(®) ~ 


$ b (0) 2 e 7ria: ' 

1 


when 5i(x) 2> 0, 
when 5i(x) <C 0. 


(A.8) 


The quantum dilogarithm is a quasi-periodic function. Explicitly, it satisfies the equations 


‘hbfr + Cb + ib) _ 1 

<f>b(T + Cb) 1 — qe 2nbx 
< hb(® + Cb + ib -1 ) _ 1 

$bO + c b ) ~~ 1 - q-ie 2 * b ~ lx 


(A.9a) 


(A.9b) 


When b is small, we can use the folllowing asymptotic expansion (see for example [76, 77]), 


log<f>[ 


(sb) 


I] (2^ib- 


k =0 


. '2k—1 B 2k ( 1/2) 
1 (2k)! 


Li2—2 fc(—e x ), 


(A.10) 


where B 2k (z) is the Bernoulli polynomial. 


B The A-period at the conifold point 


In this short Appendix we compute the A-period II^ ; (z) at the conifold point z = z c , for arbitrary 
values of m (or, equivalently, of £). The starting point of this calculation is the integral 


Z = 


1 

(27ri) 2 


/S 1 xS 1 


log P(z, w) 


d z dw 

z w ’ 


(B.l) 


where 

P(z, w) = 2(x 2 + y 2 ) + x 2 (z + z _1 ) + 2 y 2 (w + w -1 ). (B.2) 

Note that P(z,w ) is essentially the polynomial defining the mirror curve to local Fo, and the 
integral X is the logarithmic Mahler measure of this polynomial. Let us define 

4 4 

- x - y rp 

Zl 4(x 2 + y 2 ) 2 ’ Z2 4(.t 2 + y 2 ) 2 ’ 1 j 

By expanding log P(z,w) in power series in z\ t2 , we find 

1 = log (2(z 2 +! , 2 )) - £ X) r(l + ^r(i + 1)2 44. (B.4) 

n — 1 2 /c 2 / — 11 / 

If we identify the variables z\- 2 with the moduli of local Fo, we have that 



Z\ L 

m = — = 

(-) • 

(B.5) 



\yj 

and we finally obtain 

1 = 2 log y — 

l^\z c ), 

(B.6) 

where z c is given in (3.49). 
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On the other hand, the integral (B.l) was explicitly computed by Kasteleyn in section 3 of 
[78], in the analysis of the dimer model on the bipartite square lattice on the torus (see [79] for 
a nice summary of the subject.) By writing it as 

1 Z* 71 " 

1 = —g / / log (4(x 2 cos 2 u + y 2 cos 2 u/)) dwdu/, (B.7) 

n Jo Jo 

one can first perform the integral w.r.t. u/, compute the remaining integral as a power series in 
x/y = e^, and resum the resulting series in terms of the dilogarithm. One finds, 

1 = 2 log y + —Im ^ie^ ^ . (B.8) 

By comparing this to (B.6), we conclude that 

n ^ r) (^c) = ( Li 2 ( ie€ )) , (B.9) 

which is precisely what we obtained from (3.64). 
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